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ABSTRACT 


The  manager  of  a  warehouse  and  repair  facility  must  decide 
how  many  spares  of  each  of  many  different  part  types  to 
store  in  inventory  so  as  to  insure,  in  some  sense,  the 
efficiency  of  the  facility's  operations.  ^Customers  will 
bring  to  the  repair  facility  failed  parts.  If  a  spare  of 
that  part  type  is  available  from  the  warehouse,  the  customer 
will  receive  one  and  leave.  Otherwise  he  must  wait  until 
the  repair  facility  produces  a  replacement.  The  failed 
part  which  was  brought  to  the  repair  facility  will  be 
repaired  in  order  to  replenish  the  inventory  in  the  warehouse. 

Having  a  budget  constraint,  and  using  information  on  the 
failure  rates  of  the  parts,  the  relative  importance  of  each 
part  to  the  other  parts,  and  the  amount  of  time  necessary 
to  repair  a  failed  part,  the  manager  must  decide  how  to 
optimally  allocate  his  resources  in  purchasing  these  spare 
parts.  Furthermore,  he  might  also  consider  allocating  some 
of  his  funds  in  an  effort  to  increase  the  capabilities  of 
his  repair  facility. 

The  paper  first  formulates  various  mathematical  models  and 
then,  in  each  case,  shows  how  one  can  obtain  a  sequence  of 
undominated  solutions.  (  - 


INTRODUCTION  AND  SUMMARY 


Let  us  suppose  that  there  are  many  parts  of  different  types  in  use 
and  let  us  assume  that  each  of  these  part  types  is  subject  to  stochastic 
failure.  We  wish  to  set  up  a  service  facility  which  will  repair  failed 
parts  and  will  also  serve  as  a  warehouse  to  store  an  inventory  of  spare 
parts.  Let  us  say  there  are  r  different  part  types  and  that  we  have  an 
initial  inventory  of  spares  of  type  i  .  Assuming  for  the  moment  that 

k^  >  0  ,  the  following  will  occur.  Whenever  the  first  failure  of  a  part 
of  type  i  occurs,  a  spare  from  the  inventory  will  be  issued  to  the 
customer  bringing  it  in,  and  the  failed  part  will  go  to  the  repair  facility. 

Some  time  later,  another  customer  might  come  with  a  failed  part  of  type  i  . 

He,  too,  will  give  his  failed  part  to  the  repair  facility.  If  there  is  a 
spare  part  of  type  i  in  the  inventory,  he  will  immediately  receive  it. 

If  not,  he  will  have  to  wait  and  he  will  receive  the  first  part  of  type  i 
to  come  out  of  the  repair  facility.  If  a  part  of  type  i  comes  out  of 
the  repair  facility  and  no  customer  is  queued  up  waiting  for  a  part  of  this 
type,  then  the  part  goes  into  the  inventory  of  type  i  .  * 

Regardless  of  the  situation  with  respect  to  other  part  types,  if  there 
are  shortages  in  the  inventory  of  part  type  i  ,  each  customer  who  Is  waiting 
for  part  type  i — called  "customers  of  type  i" — will  receive  his  replacement 
part  on  a  first- in-first-out  (FIFO)  basis.  However,  if  a  customer  of  type 
j  is  also  waiting  for  a  part,  he  may  or  may  not  be  given  his  part  before 
those  of  type  i  receive  their  parts.  In  other  words,  there  is  no  queueing 
discipline  between  types;  only  within  types. 

As  an  example  of  a  "real-life"  system  that  this  paper  relates  to, 
consider  a  company  which  maintains  a  large  fleet  of  aircraft.  These  aircraft 
have  many  different  parts  that  are  subject  to  wear  and  even  failure.  Such  parts 
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might  be  propellers,  engines,  navigational  equipment,  flight  instruments, 
etc.  When  a  part  fails,  the  customer  brings  the  failed  part  to  the  repair 
facility  and  it  is  hoped  that  there  will  be  a  spare  part  of  the  same  type 
to  give  him  so  as  to  minimize  the  down-time  of  the  aircraft. 

The  manager  of  the  aircraft  maintenance  facility  must  determine  somehow 
the  number  of  spare  parts  of  each  type  to  store  in  his  warehouse. 

The  problem  we  wish  to  investigate  then  is  as  follows: 

How  can  one  best  allocate  his  resources  in  the  purchase  of  spare  parts 
and/or  additional  service  channels  so  as  to  minimize  the  expected  number  of 
shortages  in  the  inventory  at  some  future  time,  long  after  the  system  has 
been  in  operation? 

Chapter  1  directs  its  attention  to  the  case  in  which  only  spare  parts 
and  no  additional  service  channels  may  be  purchased.  The  method  of  approach 
in  this  chapter  is  as  follows: 

Section  1.1: 

Under  the  assumptions  given  in  Section  1.1.1,  an  expression  is  found 
for  the  probability  that  there  are  n  parts  of  type  i  in  the  system. 

By  "in  the  system,"  we  mean  in  the  repair  facility,  whether  actually  being 
worked  on  or  currently  queued  up  waiting  for  servicing. 

We  then  derive  an  expression  for  the  expected  number  of  shortages  of 
type  i  given  that  there  were  k^  parts  of  type  i  initially  in  the  inven¬ 
tory.  We  will  denote  this  by  L^(k^)  . 

Then  we  show  (under  more  general  assumptions  than  those  of  Section  1.1.1) 
that  L^fk^)  is  convex  in  k^  . 

This  is  followed  by  a  discussion  of  the  notion  of  dominance  and  the 
mathematical  optimization  problem  which  is  formulated. 

A  presentation  of  the  algorithms  of  Proschan  and  Kettelle  is  given. 

r 
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Finally,  a  generalization  is  made  to  the  case  in  which  for  any  part 
type,  there  is  only  a  subset  of  the  service  channels  which  can  work  on  it. 
The  possible  assignment  of  parts  to  service  channels  in  defined  for  a 
specific  case  and  the  methods  of  Proschan  and  Kettelle  are  shown  to  aDply. 

Section  1.2; 

An  analysis  is  presented  for  the  case  in  which  there  is  a  lag  time 
in  the  delivery  of  a  failed  part  to  the  repair  facility. 

Section  1.3: 

In  this  section  we  discuss  numerical  solutions  to  the  problem  as 
previously  defined,  but  for  systems  with  generax  arrivals  to  the  repair 
facility. 

Chapter  2  handles  the  case  in  which  funds  may  be  allocated  to  purchase 
both  spare  parts  and  additional  service  channels. 
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CHAPTER  1 


OPTIMUM  ALLOCATION  OF  RESOURCES  IN  THE  PURCHASE  OF  SPARE  PARTS 


1.1  M/M/c  Case 

1.1.1  Assumptions 

(1)  The  arrival  of  failed  parts  of  type  i  form  a  Poisson  process  at  rate 
X^  .  There  is  no  lag  between  the  time  the  part  fails  and  the  time 

it  is  brought  Into  the  repair  facility  by  the  customer.  Therefore,  the 

arrl'.al  of  all  part  types  to  the  repair  facility  forms  a  Poisson 

process  at  rate  X  «  £  X.  . 

i-1  1 

(2)  The  repair  time  for  all  parts  of  all  types  is  distributed  exponentially 

at  rate  p  ,  and  is  independent  of  the  repair  time  of  any  other  part. 

(3)  There  are  c  service  channels,  each  of  which  can  work  on  parts  of 

all  types. 

(4)  There  is  no  limit  imposed  on  the  number  of  failed  parts  that  can  queue 
up  in  the  repair  facility. 


1.1.2  Derivation  of  the  Steady  State  Probabilities 
r 

Let  X  *  l  X.  ;  c  -  number  of  service  channels;  and  p  ■  X/cp  . 
i-1 

Assume  p  <  1  .  So  far  as  the  repair  facility  is  concerned,  there  is  no 
distinction  among  parts.  Thus,  we  have  an  M/M/c  queueing  system  in  which 


M/M/c — Standard  notation  in  queueing  theory.  The  first  "M"  denotes 
exponential  interarrival  times  to  the  system.  The  second  "M"  denotes 
exponential  service  in  the  system.  The  "c"  means  there  are  c  service 
channels  in  the  repair  facility.  When  the  first  "M"  is  replaced  by 
"Cl,"  it  denotes  a  general  distribution  for  the  interarrival  times,  and 
when  the  second  "M"  is  replaced  by  a  "S,"  it  denotes  a  general 
distribution  for  the  service  time.  A  "D"  in  place  of  the  second  "M" 
denotes  deterministic  service  time. 
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the  steady  state  probability  distribution  is  given  by 


o  n! 


cc  n 
Poc!  p 


for  n  <  c 


for  n  >  c 


where  P-1  -  l  .  Let 

°  -In  nf  c! (1  -  P) 


«  number  of  type  i  in  the  system 


N  ■  number  of  all  types  together  in  the  system,  i.e., 


P(N4  -  n  |  N  -  m)  - 


:m\  -  r 


for  n  <  m 


for  n  >  m 


nh  ■ n)  ■  i  {lir)  i1-^)  p„ 


Assume  first  that  n  <  c  . 
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.uiM#  *•»  *  ~ 


Let  Oj  =  — -  ■  -i  ,  f  measure  of  the  traffic  intensity  of  all  types  other 
than  type  i  ,  and  let  ■  “ 


P(N  ®  n) 


K  I c-n  (co  )u  c-n 

P  l  - f—  +  £~r~ 

°  n!  \u-0  u!  Cl 


m-n 


l  m! 
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Now, 


l  m! 

m=c+l 


m-n 


(m  -  n)  ! 


I  m! 

m=n 


m-n 


(m  -  n) ! 
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m=>n 


_  m-n 
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(m  -  n) ! 
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Lemma  1.1: 


n ! 


(1  "  a) 


n+1 


'  "'iiwn i>it(|fc|  i  ,, u. 
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Proof : 


By  induction  on  n  .  Let  n  ■  1 


do  \1  -  0/ 


(1  -  o) 


zJlSrA 

v  2 


(1  -  o)‘ 


Now  assume  true  for  n  and  prove  for  n  +  1  ;i.e.,  assume 


(1  -  o)1 


,  n  \1  -  0  / 


/g1*1  \  j  j>  /o'*1  VI 

\1  -  o/  do  I  don  \1  -  o/J  • 


Using  the  general  formula 


[f(*)g(x)]  -  [  (jf(k)(x)g(n*k)(x) 

/fv  VcO  '  f 


where  f(x)  ■  x  and  g(x)  ■  Y~-~x  We  See  t'iat 


By  the  induction  hypothesis,  and  using  the  fact  that 


f  f f  (x)  ]  -  ~  ~~~T  [f  (x) ]  +  K  ,  we  have 
J  dxn  dx""1 


(1  -  oY 


n  7 


(1  -  0)’ 


5  ■  £[ 


- 2a!_  +  _n! - 

a  -  o)n+1  (i  -  o)n 
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d 

r  °  i  1  *■  °  il 

,  _d_ 

1 

(n  +  1) ! 

do 

L(1  -  o)  (1  -  o)  J 

■  n!  ~ 
do 

,,  Nn+1 

L(1  -  a)  J 

.n+2 
(1  -  o) 

,n+l  /  n+1 


da 


n+1 


(  on+1  \  (n  +  1)  ! 

\  l  _  0  f  *=  and  the  lemma  is  proved . 


(1  -  o)‘ 


Therefore,  for  n  <  c  , 


P(Ni  -  ">  ■  Po  ^ 


n 

p  ^ 
o  n! 


c-n  (co 

l 

u-0 


(co  )u  c-n  /  c  o  \ 

^  a. 


c-n  (coJ) 
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:o.j  c-n  /  , 

-J_  + c — / - aL 
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Now  if  n  >  c  ,  the  reader  can  see  that 


n  m-n 

P,  c-n  “  o.  c-n  . 

P(N.  =  n)  =  P  -7  ~ ~  £  m!  - - rr  -  F  p”  - - - —  . 

1  o  n!  c!  (m-n)!  o  i  c!  (1  _  }n+l 


Thus,  we  have  proved  the  following  theorem: 


Theorem  1,1: 

The  steady  state  probabilities  of  the  system  under  the  assumptions  of 
Section  1.1.1  are  given  by 


| 

t 

A 


I 


p?(" c_n  /  (C°JU  c-n  oV\  c-n  ,  1 


(1)  P(N. 


Pi  cC_n 
?°  ^  (1  -  o1)n+1 


L.1.3  The  Expected  Number  of  Shortages 


for  n  <  c 


for  n  >  c 


L.(k.)  -  l  nP(Nt  -  n  +  k  )  . 
nBl  1 


Assume  k^  <  c  -  1 


,  f  n+k.  r  . 

C_ki  0.  1  c“""kl  ( (co  )“  /-"-ki  o“\ 
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Now  the  second  term  is  equal  to 
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n+k 

cc  1  :  r  P1  1 

c!  1  “  °i  n=l  nLc(1  "  °i)-l 


k,+l 

PocC6i  y  n  0n-l 

C.d  -  0.)  jj  i 


where  0.  *= 


i  c (1  -  oi) 


.  But 


2  ‘ 


k.+l 

PocC011 

Thus,  the  second  term  equals  - r -  .  Now  If 

c!(i  -  e1)  (1  -  0l) 


ki  -  C  ’  Li^ki^ 


00 

I  n 

n**l 


n+k 


i  c-n-k. 


o  c!  n+k.+l 

(1  -  ot)  1 


But  this  is  simply  the  second  term  that  was  just  computed.  We  can  now 
write  out  a  general  formula  for  L^(k^)  • 


1  ,  n+k.  .  , 

C"ki  P.  1  C-”-ki/(co,)u  cC-"-kl  p“\ 

*  n 7r+irjT  i  l_7rr'-~Tr'<u  +  n  +  ki)l^'/ 

n=l  i  u=0 


k.+l 

«V 


c! (1  -  01)^(1  -  oi) 


for  <  c  -  1 


(2)  L^kj) 


k.+l 

«V 


0  c! (1  -  01)2(1  -  oi) 


for  k^  >  c  -  1 


•  -  *  -(tu 
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1.1. A  Proof  of  the  Convexity  of  L(k)  in  k 

/ 

Theorem  1.2: 

Under  completely  general  assumptions  on  the  distribution  of  the 
Interarrival  times  of  parts  to  the  repair  facility  and  on  the  distribution  of 
the  service  times,  and  allowing  any  number  of  service  channels  from  1  to  »  , 
the  expected  number  of  shortages  of  a  part  type  is  convex  in  the  number 
of  spares  initially  in  the  inventory. 

Proof ; 

Let 

k  -  0  of  spares  initially  in  the  inventory 

P  *  P{n  parts  in  the  system  (steady  or  transient  state)} 
n 


L(k) 


00 

l 

n*l 


nP 


n+k 


A  L(k)  -  L(k  +  1)  -  L(k) 

A2L(k)  -  L(k  +  2)  -  2L(k  +  1)  +  L(k) 


"  ^  ntPn+k+2  "  2Pn+k+l  +  Pn+k‘ 

n=l 
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Note  that  along  the  indicated  diagonals,  the  sum  is  zero.  Thus,  the 


total  sum  is  merely  >  0  .  // 


1.1.5  The  Optimization  Problem 


..1.5.1  Notion  of  Dominance 


For  each  allocation  Jt  *=  (k^,k2»  ...»  kf)  ,  there  is  a  value  (in  the 
sense  of  "worth")  and  a  cost  associated  with  it,  denoted  by  V(k)  and 
C(k)  respectively. 

Definition:  Allocation  k.^  is  said  to  dominate  allocation  k^  iff 
V(kx)  >  V(k2)  and  Cflcj)  <C(k2)  . 


Definition:  Allocation  Jtj  is  sa id  to  strictly  dominate  allocation  k^ 
iff  at  least  one  of  the  above  inequalities  is  strict. 

Definition;  Allocation  k^  is  said  to  be  undominated  if  there  does  not 
exist  an  allocation  k^  which  strictly  dominates  it. 

Definition:  A  sequence  S  of  undominated  allocations  is  called  a  complete 
fctnily  of  undominated  allocations  iff  for  every  allocation  It  , 
either  k  c  S  or  there  exists  an  element  of  S  which  dominates 


Definition:  A  sequence  S  of  undorainated  allocation,  which  is  not  a 

complete  family,  is  called  an  incomplete  family  of  undominated 
allocations. 


1.1. 5. 2  Mathematical  Formulation  of  the  Optimization  Problem 


Let  us  say  that  each  part  of  type  i  has  a  cost  c^  .  Thus,  allocation 

r 


It  has  a  cost  associated  with  it  of  C (Jk )  *  £  c.k.  . 

i«l 


Furthermore,  let  us 


■  ' 
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say  that  there  is  a  constant  v^  associated  with  parts  of  type  i  which 
is  a  measure  of  the  relative  importance  of  part  type  i  which  is  a  measure 
the  relative  importance  of  part  type  i  to  the  other  parts  ,  For  examnle, 
v  might  be  the  cost  in  dollars  for  a  shortage  of  part  type  i  . 

We  wish  to  construct  a  sequence  of  undominated  allocations  to  the 
problem  with  V(k)  ■  -W(k)  ,  where 

r 

(3)  W(k)  -  l  v  L  (k.) 

i-1  11  1 

and 

r 

(4)  CQc)  «  l  c  kt  . 

i»l 

Having  this  sequence  of  undominated  allocations,  one  knows  that  if 

* 

.k  is  an  element  of  the  sequence  then  for  any  allocation  k  with 
W (jc * )  <  W(k  )  ,  then  it  must  be  the  case  that  C(k')  >  C(k  )  , 

If  the  sequence  is  an  incomplete  family  of  undominated  allocations, 
the  following  problem  arises.  If  exactly  b  dollars  are  allotted  to  be 
spent  purchasing  spare  parts,  and  if  in  constructing  the  incomplete  family 
of  undominated  allocations  one  finds  no  solution  with  associated  cost 
exactly  equal  to  b  ,  but  finds  the  two  solutions  k^  and  k^  such  that 
CQcj)  <  b  <  CQcj,)  .  and  C(]^)  is  closest  to  b  from  below  while  CQ^) 
is  closest  to  b  from  above,  then  the  optimal  allocation  vector  to  the  problem 

Minimize  W(Jc) 

(5) 

Subject  to  C(k)  <  b 

t 

might  be  k^  ,  but  it  is  possible  that  there  is  an  undominated  allocation 
which  the  algorithm  skipped  over  and  this  allocation  has  ccrt  between  CQc^) 
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and  b  .  In  any  case,  if  we  choose  allocation  as  the  optimum  vector 

of  allocations  to  solve  (5)  above,  then  our  error  in  the  functional  value 
is  no  greater  than  WQcj)  -  W(k^)  and  the  error  in  the  cost  is  no  greater 
than  b  -  C(k^)  . 

If,  however,  the  generated  sequence  is  a  complete  family  of  undomlnated 
allocations,  the  problem  just  stated  does  not  arise.  That  is,  the  optimal 
vector  of  allocations  to  (5)  is  . 

1.1.6  Two  Algorithms  to  Generate  Sequences  of  Undominated  Allocations 
1.1.6. 1  Proschan’s  Algorithm 

The  following  algorithm  is  taken  directly  from  reference  [1].  The 
notation  and  wording  has  been  only  slightly  changed  to  be  more  meaningful 
to  this  problem. 

Start  with  the  cheapest  cost  allocation  k  *  (0,0 . 0)  . 

Obtain  successively  more  expensive  allocations  as  follows.  If 
our  present  allocation  is  k  ,  we  determine  the  index,  say  i  , 
for  which 


^ 4  VV 


is  minimum  over  i  ■  1,2,  ...,  r  .  (If  the  minimum  is  achieved  for 
more  than  one  value  of  the  index,  choose  the  lowest  among  these.) 

Then  the  next  allocation  is  /k^,k2>  ...»  k.  _^,k  +l»k^  +^,  ...,  k  \  ; 

'  o  o  o  ' 

that  is,  we  have  added  a  single  unit  of  the  io-th  type  to  k  . 

The  "A"  sign  in  the  above  expression  is  defined  to  be 
AL^)  -  +  1)  -  L^)  . 


It  can  be  shown  that  for  the  function  L^(k^)  as  given  in  Equation  (2) 
(where  the  subscript  "i"  has  been  dropped), 
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1c— k— 1  n+k  c-n-k-1  f  /  v  .  .u 

Jx  "  ^  J0  [M-  -  l)  ^ 


ec9lt*1 


c!  (l  -  e)(i  -  o) 


for  0  <  k  <  c  -  2 


(6)  AL(k) 


r  ccek+1 
0 


c!  (1  -  9) (1  -  o) 


for  k  >  c  -  1  , 


r 

Proschan  has  shown  that  if  W(k)  ■  £  *8  convex  in  k  » 

i*l 

then  the  previously  stated  algorithm  produces  an  Incomplete  family 
of  undominated  allocations.  Since  Section  1.1.4  shows  that  L^k^)  is 
convex  in  k^  ,  we  know  that  W(k)  is  convex  in  k^  .  Thus,  the  algorithm 
applies. 

It  should  be  noted  that  Proschan's  algorithm  will  handle  multiple 
constraints.  But  since  this  does  not  add  to  the  purpose  of  this  paper, 
a  discussion  of  multiple  constraints  is  omitted. 


1.1. 6. 2  Kettelle's  Algorithm 

In  this  section,  we  construct  a  complete  family  of  undominated  allo¬ 
cations  for  the  problem  in  which  W(k)  and  C(k)  are  given  by  (3)  and 
(4),  respectively.  The  method  employed  is  that  of  Kettelle. 

The  algorithm  works  progressively  in  that  it  first  obtains  a  complete 
family  of  undominated  allocations  for  the  subsystem  in  which  only  two  of 
the  r  part  types  are  present.  Then  it  takes  another  two  of  the  r  part 
types  and  constructs  a  complete  family  of  undomlnated  allocations  for  this 
subsystem.  Now  it  takes  these  two  complete  families  and  joins  them  in 
such  a  way  as  to  generate  from  them  a  complete  family  for  the  subsystem 
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containing  all  four  of  the  part  types.  The  method  proceeds  In  that  manner 
so  that  at  each  stage,  a  complete  family  is  generated  for  successively 
larger  subsystems  from  two  smaller  subsystems.  There  is  no  requirement 
that  the  subsystems  be  of  equal  size.  Nor  is  there  any  requirement  of 
convexity  of  W(k)  as  there  was  in  Proschan's  algorithm. 

For  the  exact  algorithm  at  each  stage,  the  reader  is  referred  to 
Kettelle  [3]. 

It  should  be  noted  that  Kettelle' s  algorithm  is  applicable  only  when 
there  is  one  cost  constraint.  Proschan  and  Bray  [7]  have  extended  the 
method,  however,  to  multiple  constraints. 

As  in  the  previous  section,  a  discussion  of  the  multiple  constraint 
case  is  omitted. 


1.1.7  A  Simple  Extension 

Let  ya  represent  the  i-th  service  channel. 

Let  T  *  {all  service  channels}  -  { |  i  ■  1,2 . c) 


T(i) 


service  channel  j  can  service  part  type  i} 


Assume  the  following: 

(1)  T(i)  i  0  for  all  i  -  1,2,  ....  r  . 

(2)  For  all  J  ,  there  exists  i  such  that  Yj  c  T(i)  . 

(3)  T(i)  c  t  and  T(i)  +  T  for  all  i  -  1,2,  ....  r  . 

(4)  If  T(i)  0  T(j)  t  0  ,  then  T(i)  -  T(j)  for  all  i,j  -  1,2 r  . 


These  assumptions  result  in  a  unique  partition  of  T  into  q  disjoint 


subsets, 
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T  *=  T,  +  T-  +  ...  +  T  . 

12  q 

This  partition  then  separates  the  entire  system  into  q  completely 
Independent  queueing  systems. 

Let  each  subset  have  a  parameter  associated  with  it  which 

is  the  rate  at  which  all  Yj  €  produce  exponential  service. 

We  now  wish  to  find  a  sequence  of  undominated  allocations  to  this 
problem  with 


W(k)  -  l  l  v  L  (k  ) 

i-1  jcT1  J  J  J 

(7) 

C(k)  -  ?  I  c  k  . 

i-1  jcT1  J  J 


It  is  easy  to  see  that  the  form  of  (7)  in  no  way  complicates  the 
method  of  solution  as  outlined  in  either  algorithm  of  Section  1.1.6. 

Thus,  the  methods  of  both  Proschan  and  Kettelle  readily  apply  to 
this  extension  of  the  problem. 

1.2  Lag  in  Delivery 

Suppose  that  upon  failure  of  a  part  type,  the  failed  part  does  not 
apoear  Instantaneously  at  the  repair  facility.  Rather,  let  us  assume 
that  there  is  a  lag  time  D  ,  a  random  variable  with  general  distribution 
H(x)  (assumed  different  for  each  part  type). 

As  in  Section  1.1,  let  us  assume  that  the  occurrence  stream  of 
failures  of  part  type  i  forms  a  Poisson  process  at  rate  .  We  are 
Interested  in  determining  the  arrival  stream  to  the  repair  facility. 
Dropping  the  subscript  i  for  the  moment,  let  us  define 

Sj  ■  time  the  j-th  failure  (of  part  type  i)  occurs 
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Dj  =  lag  time  for  j-th  failure  before  it  appears  at  the  repair  facility 

Yj  *  time  of  the  j-th  arrival  to  the  repair  facility  (not  necessarily  the 
same  as  the  time  the  j-th  failure  arrives  at  the  repair  facility). 

The  following  figure  represents  a  possible  outcome  in  the  failure 
and  arrival  pattern. 


FIGURE  1.1:  POSSIBLE  OUTCOME  IN  THE  FAILURE 
AND  ARRIVAL  PATTERN 

We  must  determine  the  distribution  of  Y^+1  -  Y^  ,  the  interarrival 
time  of  failed  parts  to  the  repair  facility.  To  do  this,  consider  an 
imaginary  M/G/09  queueing  system.  The  occurrence  of  a  part  falling  will 
be  an  arrival  to  the  imaginary  system.  In  this  imaginary  system,  the 
part  goes  immediately  into  service,  no  matter  how  many  parts  have  already 
arrived  (i.e.,  failed);  thus,  there  are  an  infinite  number  of  servers  in 
the  imaginary  system.  The  service  is  analogous  to  the  lag  time  D  .  When 
the  part  leaves  the  imaginary  system,  it  has  completed  its  lag  time  and 
thus  it  arrives  at  our  actual  repair  facility. 

Therefore,  we  can  say  that  the  arrival  stream  to  our  repair  facility 
is  the  output  stream  of  an  M/G/00  queueing  system,  where  the  service 
distribution  is  H(x)  . 

Mirasol  [4]  has  shown  that  the  output  stream  of  an  M/G/09  queueing 
system  is  a  nonhomogeneous  Poisson  process  such  that  if  iKt,T)  *  number 
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of  customers  to  leave  during  (t,t  +  T]  ,  then 


t+T 


i|)(t,T)  -  Poisson  ‘  /  G(x)dx I  . 


i 

i 

Taking  the  limit  as  t  -*■  ®  of  the  parameter  shows  that  in  steady 
state,  the  output  stream  is  a  Poisson  process  at  rate  X  ,  the  input  rate. 

Thus,  even  under  these  more  general  conditions,  the  interarrival 
times  of  failed  parts  to  the  system  is  exponential  at  rate  X  .  Introducing 
the  subscript  i  again,  we  can  say  that  the  arrival  stream  of  parts  of 
type  i  to  the  repair  facility  is  a  Poisson  process  at  rate  X^  . 

Let  us  say  that  for  part  type  i  ,  D  -  H^(x)  and  let  us  let 


-  J  xdH^x)  ,  the  expected  value  of  D  for  the  i-th  part  type.  Let 

w^(k^)  "  t*ie  exPect8C*  waiting  time  a  customer  incurs  between  the  time  the 
failed  part  arrives  to  the  repair  facility  and  the  time  a  replacement  is 
available  for  the  customer,  given  that  parts  of  type  i  are  initially 

in  the  inventory.  Finally,  let  +  w^Ck^)  ,  the  expected  total  time 

a  customer  of  type  1  must  wait  from  the  moment  his  part  fails  to  the 
time  when  he  has  a  replacement  from  the  repair  facility. 

Let  us  assume  the  following: 


(1)  There  are  an  infinite  number  of  service  channels. 

(2)  The  service  time  required  to  repair  a  part  of  type  i  is  a  random 
variable  with  general  distribution  G^(x)  . 

Under  these  assumptions,  Proschan  [6]  has  shown  that  the  expected 
number  of  shortages  (in  steady  state)  of  parts  of  type  i  ,  given  k^ 
of  them  are  initially  in  the  inventory  is  given  by 


I 
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(8) 


Li(kl) 


l 

J-k±+l 


(j  - 


(X1n1)^ 

_ 


where 


00 


xdG^(x)  . 


In  order  to  find  w^(k^)  ,  observe  the  following: 

We  define  an  imaginary  queueing  system  (not  the  one  previously 
defined  in  this  section)  in  which  the  steady  state  probabilities  I  are 
given  by 


I 

n 


i 

I  P 
-0 


n  ■  0 


n  >  0 


where  is  the  steady  state  probability  that  there  are  n  of  type  i 

in  the  system. 

Let  J  *  the  expected  number  in  the  imaginary  system.  Using  the 
well  known  relation  "L  -  Aw"  (see  [2])  where  L  represents  the  expected 
nutaber  in  the  system,  A  the  arrival  rate  of  customers  to  the  syscem,  and 
w  the  expected  waiting  time  a  customer  incurs  in  the  system,  we  have  that 
the  expected  waiting  time  in  our  imaginary  system  is  J/A^  ,  But 
J  ■  L^(k^)  .  Therefore, 


(9)  wi(ki)  ■=  Li(k1)/A1  . 

If  x ^  ,  as  defined  earlier  in  this  section,  is  the  expected  total 
time  a  customer  of  type  i  must  wait  from  the  moment  his  part  fails  to  the 
time  when  he  has  a  replacement  from  the  repair  facility,  and  if  the 


umrnmm 


* 
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probability  that  an  arbitrary  customer  is  of  type  i  is  X^/ X  ,  then  the 
expected  total  waiting  time  of  an  arbitrary  customer  is  given  by 


r 

I 

i-i 


Vi 


( * 

u«i 


Vi + 


r 

l 

i=l 


XiwJ (k 


‘V1  • 


Introducting  the  constants  as  was  done  previously,  and  using  (9), 


(10) 


W(k) 


Vi°i 


r 


+  l 

i-1 


ViW 


Since  be  Section  1.1.4  L^(k^)  is  convex  in  ,  it  follows  that  W(k) 
is  convex  in  Jc  .  This  being  the  case,  both  the  algorithm  of  Proschan  and 
that  of  Kettelle  apply. 

It  is  interesting  to  note  that  whether  one  applies  one  algorithm  or 
the  other,  the  resulting  sequence  of  undominated  allocations  will  be 
identical  to  the  resulting  sequence  if  he  had  applied  the  same  algorithm 
to  the  problem  of  minimizing  the  weighted  sum  of  the  expected  number  of 
shortages.  That  this  is  true  is  evident  from  (10). 

1.3  Systems  with  General  Arrivals 
1.3.1  GI/D/™  Case 

Assume  the  distribution  of  the  interarrival  times  to  the  repair  facility 
to  be  G^(x)  for  the  i-th  part  type.  Assume  further  that  the  service 
distribution  for  this  part  type  is  deterministic  at  value  d^  .  Finally, 
assume  that  there  are  an  infinite  number  of  servers  available  to  work  on 
failed  parts.  Thus,  the  underlying  queueing  system  is  assumed  to  be  GI/D/00  . 

Let  us  solve  for  the  steady  state  probabilities  of  this  system. 

Let  Pn(t)  =  parts  of  type  i  in  system  at  time  t)  .  Then, 
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P  (t  +  d  )  +  l  P. (t)*P{n  parts  of  type  i  arrive  during  the  interval 

1  j-0  2 

(t,t  +  d^]}  .  Since  the  infinite  sum  above  is  one,  we  have  as  t  00 

Pn  «=  P{n  parts  of  type  i  arrive  during  an  interval  of  length  d^,  "long 
after"  the  system  has  been  in  operation}  . 

Assume  initially  that  n  >  1  .  Let  us  condition  on  the  arrival  of  the 
first  of  these  n  arrivals.  Suppose  it  occurs  x  units  of  time  after  the 
interval  has  started,  and  since  n>l,0<x<d^.  Then  during  the 
remainder  of  the  interval,  d.  -  x  ,  exactly  n  -  1  arrivals  must  occur. 
Figure  1.2  illustrates  this  discussion. 


first  n-1  arrivals 
arrival  during  this 


interval 

FIGURE  1.2:  ARRIVALS  DURING  INTERVAL  OF  LENGTH  dL 

Now  the  probability  of  exactly  n-1  arrivals  during  the  interval  of 
length  d^  -  x  is  given  by 

o<n-i>(d1  -  x)  -  <=;">  (d± .  x> 

where  C is  the  n-fold  convolution  of  G.  with  Itself. 

In  order  to  determine  the  distribution  of  x  ,  we  note  that  x  is  the 
"residual  life"  or  "excess  time"  until  another  arrival  occurs.  It  is  known 


Mtl 


20 


(see  [8])  that  in  steady  state,  the  excess  time  has  the  "equilibrium 

0 

distribution"  of  G^  ,  denoted  by  ,  where 

tm 

C,(t)  «  -f-  /  [1  -  G  (u)]du 

i  J0 


and  p.  ■  J  udG  (u)  . 

1  0 

Thus,  unconditioning  the  above  probability  expression,  we  have  for 
n  >  1  , 


(11. a)  P 


n 


J  [Gf"1)<di  -  *>  -  G^n>(d1  -  *>]< 


(x)  . 


For  n  •=  0  , 


(11. b) 


p0  -  1  -  o'M,)  . 


The  expected  number  of  shortages  of  part  type  i  is  given  by 


Since  by  Section  1.1.4,  L^(k^)  is  convex  in  ,  we  know  that 

r 

W(k)  *  l  v.L, (k, )  is  convex  in  k  and  either  of  the  algorithms  of 
i=l  ? 

section  1.1.6  apply  to  generate  a  sequence  of  undominated  allocations. 

As  done  in  Section  1.2,  we  could  apply  "L  ■  XW"  if  we  desired  to 
minimize  the  expected  waiting  time  instead  of  the  expected  number  of  shortages. 
However,  as  was  the  case  in  Section  1.2,  the  resulting  optimal  allocation 
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vector,  jc  ,  would  be  Identical  to  the  optimal  allocation  vector  obtained 
in  minimizing  the  expected  number  of  shortages. 

It  should  be  mentioned  that  even  though  one  might  have  an  analytical 
expression  for  G^(x)  ,  he  will  not,  in  general,  be  able  to  get  a  working 
expression  for  G^n^ (x)  ,  or  for  that  matter.  Equations  (11. a)  and  (11. b). 
Therefore,  he  will  have  to  resort  to  a  numerical  solution  to  solve  (11). 
This,  though  possibly  time  consuming  even  on  a  high  speed  computer,  should 
not  present  too  great  a  difficulty. 

1.3.2  GI/G/c — Simulation  Studies 

In  some  cases,  it  will  not  be  feasible  to  make  the  simplifying 
assumptions  of  the  previous  sections  with  regards  to  the  arrival  stream 
and/or  the  distribution  of  service  times.  Unless  one  is  able  to  derive 
an  analytical  expression  for  {P^}  under  his  particular  set  of  assumptions, 
he  will  be  forced  to  use  simulation.  If  one  uses  simulation,  there  is,  of 
course,  no  restriction  on  the  number  of  service  channels  he  may  specify 
for  the  model. 

Having  obtained  the  {P^ }  through  simulation,  one  may  then  numerically 
proceed  with  the  analysis  to  obtain  an  expression  for  W(k)  ,  and,  by 
virtue  of  Section  1.1.4,  he  may  apply  either  of  the  algorithms  of  Section 


1.1.6. 
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CHAPTER  2 

OPTIMUM  ALLOCATION  OF  RESOURCES  IN  THE  PURCHASE  OF  SPARE  PARTS 
AND  ADDITIONAL  SERVICE  CHANNELS 


2.1  Introduction 

Consider  the  problem  of  Section  1.1  with  the  following  change. 

Instead  of  being  allowed  to  spend  the  resources  only  on  spnre  parts  of 
each  type,  we  are  now  allowed  to  allocate  an  unspecified  proportion  of 
the  money  to  purchase  (or  create)  additional  service  channels.  How  then 
should  we  allocate  our  resources? 

2.2  Algorithm 

We  assume  wc  start  with  c  service  channels  and  the  cost  of  these  is 
zero.  Let  s  be  the  number  of  additional  channels,  and  let  Q(-)  be  the 
cost  to  purchase  them.  Two  assumptions  on  Q(s)  are  made: 

(1)  AQ(s)  *  Q(s  +  1)  -  Q(s)  >  0  for  all  s  -  0,1,2,  ... 

(2)  lim  Q(s)  *  * 

8-**° 

In  view  of  Assumption  (2),  and  since  it  is  not  unreasonable  to  assume 
a  space  constraint,  let  us  say  that  under  no  circumstances  will  we  purchase 
more  than  q  additional  service  channels. 

We  now  wish  to  construct  a  complete  family  of  undomlnated  allocation3 
to  this  generalized  problem.  The  algorithm  that  we  will  employ  is  as 
follows: 

I.  Using  Kettelle's  algorithm,  assume  c  +  s  service  channels,  and 
construct  a  complete  family  of  undominated  allocations.  Do  this 
for  s  ■  0,1,2,  ...,  q  . 

II.  Construct  the  following  table: 
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if  of 
Service 
Channels 

Zero-th 

Allocation 

First 

Allocation 

n-th 

Allocation 

c 

<2.0  ,M00) 

^l^Ol’^W 

^0n,b0n,W0n)  *  ‘  * 

c  +  1 

<2-bio'V 

%l*bll*Wll) 

(-ln,bln,Wln;-  '  ' 

c  +  2 

• 

^’b20,VW 

• 

%l,b21’W21) 

• 

^2n,b2n’^2n^  ’  *  * 

• 

• 

• 

C  +  q 

• 

• 

<-°-bqo-V 

• 

• 

vvv 

• 

• 

(k  ,b  ,W  ).  ,  . 

^qn  qn  qn 

where  k^n  is  the  n-th  member  of  the  complete  family  of 

undominated  allocations,  assuming  c  +  i  service  channels, 

b.  is  the  cost  of  the  allocation  vector  plus  the  cost 

In  in 

Q(i)  of  i  additional  service  channels,  and  is  the 

functional  value  W(k^)  . 

III.  t  -  1 

<k,..bt,Wt)  -  (0,0, WQ0) 


IV. 


j0  -  1  J  ji  *=  °  for  all  i  -  1,2,  ... 

i  ■  0.  Compare  W  to  W, 

C  13i 

(a)  If  Wt  <  ,  set  -  j  ±  +  1 

(b)  If  W  >  W  .  ,  continue 

t  iJj 


»  q  • 


and  repeat 


IV- (a) 


i  -  i  +  1 

If  i  >  q  ,  go  to  V 
Otherwise,  go  to  IV- (a)  . 
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V.  t  -  t  +  1 

Let  s  be  the  Index  such  that 


Min  (b 


ij. 


If  there  is  a  tie,  choose  s  such  that  W  .  is  the  minimum  of 
all  those  with  equal  cost. 

Set 


j  -  j  +  1 

Js  Js 


Go  to  Step  IV 


The  resulting  sequence  represents  a  complete  family  of  undominated 
allocations  for  the  generalized  problem,  as  proved  in  Theorem  2.1. 


Theorem  2.1; 

The  allocations  obtained  using  I-V  constitute  a  complete  family  of 
undominated  allocations. 

Proof : 

First  we  will  show  that  all  allocations  found  by  the  procedure  are 

t 

undominated. 

Assume  the  contrary  .  Say  there  exists  m,t  such  that  W  <  W  and 

m  t 

<  bfc  ,  and  both  of  these  solutions  were  picked  by  the  procedure. 

Steps  IV  and  V  tell  us  this  is  not  possible  because  if  allocation  t 
were  chosen,  then  any  succeeding  choice  would  have  higher  cost  and  lower 
functional  value. 

Now  suppose  m  is  such  that  allocation  m  was  not  chosen  by  the 
procedure.  And  let  us  say  that  allocation  m  dominates  allocation  z  which 
the  procedure  chose.  Denote  this  by 
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A  >  A  . 
m  z 


By  the  definition  of  dominance  in  Section  1.1.6. 1, 


A  >  A  <=»  W  <  W  and  b  <  b 
m  z  m  z  m  z 


Let  the  sequence  of  chosen  allocations  be 


S  «  {ArA2,  ....  An,An+1,  ....  Az,  ...} 


Let  n  be  such  that  b  <  b  <  b  . 

n  m  n+1 

Since  W  <  W  ,  and  W  <  W  ,  we  have  W  <  W  , ,  which  implies  A 
z  n+1  m  z  m  n+1  m 

dominates  all  allocations  from  A  . ,  to  A  . 

n+1  z 

When  A  was  chosen  as  the  allocation  after  A  ,  either  A  was 
n+1  n  m 

"checked"  or  it  was  not.  By  "checked,"  we  mean  that  in  Step  V,  A^  was 
one  of  the  q  +  1  allocations  which  were  compared.  If  Am  was  Indeed 
checked,  then  by  virtue  of  this  fact. 


(13) 


V  <  W  and  b  >  b  . 
m  n  m  n 


Since  A  ,,  was  also  checked, 
n+J. 


(1A) 


W  , ,  <  W  and  b  , .  >  b 
n+1  n  n+1  n 


Since  by  Assumption  Affl  >  A^+^  , 


(15) 


W  <  W  <  W  and  b  <  b  <  b  . 
m  n+1  n  n  m  n+1 


But  if  this  is  true,  A  would  have  been  chosen  over  A  ,  and 

m  n+1 

A  would  never  have  been  chosen.  Thus,  A  could  not  have  been  checked. 
n+1  m 
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Now  if  A  war  uot  checked,  there  are  two  possibilities, 
m  r 


(i)  was  already  discarded  as  a  possible  choice. 


(ii)  A^  has  not  yet  been  checked. 


Suppose  (i)  is  the  case.  But  since  W  <  W  ,  A  could  not  have  been 
rr  m  n  m 

discarded  because  Step  IV  of  the  procedure  will  discard  Am  only  if 

W  <  W  Therefore,  (i)  is  not  possible. 

n  -  m 

Suppose  (ii)  occurs.  This  means  there  exists  A^  in  the  same  row  of 

the  table  as  A  such  that 
m 


(16) 


W  <  W  and  b  >  b  . 
m  v  m  v 


Le;  as  say  that  is  being  checked  as  a  possibility  to  become  An+j 

First  note  that  A  could  not  have  been  the  one  chosen  for  A  . , 
v  n+1 

because  A  >  A  , ,  ,  but  A  ^  A  . 

m  n+1  in  v 

Now  since  Ay  was  checked,  the  algorithm  says 


(17) 


W  <  W  and  b  >  b  . 
v  n  v  n 


And  since  A  . ,  was  checked  and  chosen  over  A  , 
n+1  v  ’ 


(18) 


W  <  W  and  b  .  >  b 
n+1  n  n+1  n 


and 

(19) 


b  , ,  <  b 
n+1  v 


Inequalities  (16)  -  (19)  imply 


b  >  b  >  b  b 
m  v  n+1  n 
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But  we  have  already  said  that  b  >  b  >  b 
J  n+1  m  n 

Thus,  (ii)  is  not  possible  and  we  have  shown  that  no  allocation  in 

the  table  dominates  a  chosen  allocation. 

In  order  to  prove  that  all  chosen  allocations  are  undominated,  it 

remains  to  show  that  if  we  have  an  allocation  A  which  Kettelle's 

m 

algorithm  did  not  choose  when  Step  I  was  performed,  then  could  not 

dominate  any  chosen  allocation  in  S  . 

Assume  otherwise.  That  is,  there  exists  A  e  S  such  that  A  >  A  . 

z  m  z 

Since  Kettelle's  algorithm  did  not  choose  A  ,  there  exists  A  ,  an 

m  t 

allocation  with  the  same  number  of  additional  service  channels  as  in 

allocation  A  ,  such  that  A_  ^  A 
m  t  m 

Now  Afc  i  S  ,  for  if  it  were,  then  since 


A,>  A>  A 
t  m 


z 


» 


we  would  have  two  chosen  allocations  in  S  ,  one  of  which  dominates  the 
other.  This  is  not  possible. 

But  even  if  Af  were  not  chosen,  we  have  just  proved  that  no  allocation 
in  the  table  (such  as  A^.)  dominates  a  chosen  allocation  in  S  (such  as 

V* 

Now  we  must  show  the  jonverse;  i.e.,  if  A  i  S  ,  there  exists 

'  ’  m 

A  c  S  such  that  A  V  A 
z  z  m 

Since  A  i  S  ,  there  had  to  be  some  point  at  which  it  was  eliminated 
tn 

from  further  consideration  to  become  a  member  of  S  .  Let  us  sa”  that  A 

;  m 

was  eliminated  when  we  were  looking  for  Av+^  ,  having  just  chosen  Ay  . 

If  A  was  eliminated,  then  Step  IV  says  that  W  <  W  and  since  A 
m  v  -  m  tn 

was  checked ,  b  >  b 
m  v 

Thus,  A  >  A  and  there  exists  a  chosen  allocation  which  dominates 
’  v  m 


A  .// 
m 


* 

K 

\ 


■**»&*"  nmmmm tt*+v*u* 
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It  should  be  noted  that  the  chosen  sequence  of  allocations  might 
have  elements  which  would  be  dominated  by  an  allocation  in  which  the  number 
of  service  channels  exceeds  c  +  q  .  But  since  we  have  said  that  c  +  q 
is  the  maximum  number  of  service  channels  which  wp  will  consider,  there  is 
no  need  to  be  concerned  over  this. 
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